Measures of globalization based on cross-correlations of world financial indices. 
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The cross-correlation matrix of daily returns of stock market indices in a diverse set of 37 countries 
worldwide was analyzed. Comparison of the spectrum of this matrix with predictions of random 
matrix theory provides an empirical evidence of strong interactions between individual economies, as 
manifested by three largest eigenvalues and the corresponding set of stable, non-random eigenvectors. 
The observed correlation structure is robust with respect to changes in the time horizon of returns 
ranging from 1 to 10 trading days, and to replacing individual returns with just their signs. This last 
observation confirms that it is mostly correlations in signs and not absolute values of fluctuations, 
which are responsible for the observed effect. Correlations between different trading days seem to 
persist for up to 3 days before decaying to the level of the background noise. 

PACS numbers: 89.65.Gh, 89.75.Fb, O5.40.Ca, 89.7Q.-Fc 



In spite of the tremendous importance that current 
public opinion places on issues of globalization of the 
world's economy, its sources and consequences remain 
poorly understood. Large downturns and collapses of 
the economic and financial situation in one country are 
routinely blamed on recent events in other countries. 
This point of view is reinforced by sensational newspaper 
headlines like "Latin American markets catch the asian 
flu" . The level of globalization of a diverse set of 50 de- 
veloped countries and key emerging markets worldwide 
was recently measured and reported in the A.T. Kear- 
ney/Foreign Policy Magazine 

Globahzation IndexTM [|. 
The factors selected to contribute to this index are ex- 
tremely diverse and include among other things volumes 
of inward- and outward-directed foreign investments, the 
amount of international travel and phone calls, number 
of servers of the World Wide Web, etc. Among other 
things globalization is expected to manifest itself in the 
dynamics of financial indices of stock markets in differ- 
ent countries. Indeed, it is reasonable to expect that a 
significant coupling of the economy of a given country to 
the rest of the world (e.g through foreign investments), 
would make its stock index more susceptible to changes 
in the world economic climate. 

In this work we suggest a simple measure of the level 
of financial globalization of a given country based on 
the analysis of cross-correlations between stock market 
indices in different countries and regions of the world. 
The main object of our study is the N x N empirical 
correlation matrix Cy of index price fluctuations in a 
large number of individual countries {N — 37 in our 
study). The matrix is constructed by applying the for- 
mula Cij = j^^t^i X to the set of normal- 
ized local currency returns of individual indices, recorded 
over a period of T trading days. A return SXi{t) of the 
stock index Si{t) with the time horizon At is usually 
defined as SX-{t) = lIiS^{t + At) - \IiS^{t) ~ {S^{t + 
At) — Si{t)) / Si{t) . Different markets are characterized 
by different volatilities of their stock market indices. In 



order to be able to detect similarity in the pattern of 
returns in different countries on needs to exclude volatil- 
ity effects by using normalized returns Axi{t). Axi{t) 
is constructed by offsetting each dXi{t) by its empirical 
average value (dXi), and normalizing it by its empirical 
variance (volatility): Sxi(t) = f^'^!^"^,*^^'^, . The ma- 

y/{SXf)-{SX,)^ 

trix Cij defined using Sxi (t) has the property that in the 
absence of correlations and in the limit T it is just 
the unity matrix. However in real life the number of trad- 
ing days T in one's dataset is always finite. As a result 
an empirically measured matrix dj is always dressed by 
a substantial amount of noise. It is exactly the task of 
separating any real correlations present in the signal from 
this spurious noise, that makes the analysis of real world 
data highly non-trivial. Mainstream economics literature 
was mostly devoted to a detailed analysis of these corre- 
lations for just a pair of stock indices, e.g. those of New 
York and Tokyo Stock Exchanges ||^ , but even if a large 
correlation matrix was considered usually little effort was 
made to reliably separate the signal from the noise. As 
is common in statistics, the job of looking for correlation 
patterns among several noisy signals gets much simpler 
when the number of signals is large. In our data this cor- 
responds to a large number of country indices N. The 
spectral analysis of the correlation matrix followed by a 
comparison of the spectrum with predictions of the Ran- 
dom Matrix Theory is a useful tool which allows one to 
detect even weak correlations between multiple signals. 
This method was recently successfully applied towards 
finding reproducible correlations between price fluctua- 
tions of hundreds of stocks traded on the stock exchange 
of a single country ||], and two countries Q. In this 
work we go one step further and apply the techniques pi- 
oneered in 1^ to a large number of stock market indices in 
a geographically diverse set of countries. An alternative 
method of analysis of the financial data from the matrix 
of correlation coefficients by constructing the Minimal 
Spanning Tree (MST) was described in and recently 
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applied to cross-correlations of world financial indices in 
1^. Authors of this work also found a non-trivial corre- 
lation pattern manifested by a strong regional grouping 
of indices in the MST. 

The raw data we had at our disposal consists of the 
daily open, high, low, and close prices of leading mar- 
ket indices (one per country) in 15 European, 14 Asian, 
and 8 North and South American countries. We first cal- 
culated the daily open-to-close returns of each of these 
indices. We further selected from our set only those trad- 
ing days for which we had a valid record for each and 
every country on our list. All data have gaps in them 
e.g. due to national holidays, when a particular mar- 
ket was closed. That left us with precisely 226 trad- 
ing days approximately uniformly distributed between 
April 28, 1998 and December 20, 2000. Each of these 
remaining daily returns was normalized in such a way 
that = 0' and = 1- The his- 

togram of all 37 eigenvalues of the correlation matrix 
Cij, shown in Fig. 1(a), revealed that the majority of 
eigenvalues are consistent with a null hypothesis of inde- 
pendent identically distributed Gaussian variables Sxi (t) . 
The prediction for the eigenvalue density 
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given for this situation by Random Matrix Theory |^] 
reasonably agrees with our data below A ~ 1.3. This 
formula, derived in the limit of very large T and N , pre- 
dicts sharp lower and upper cutoffs, A± = a'^{l + [N/T)± 
2y^N/T), in the eigenvalue density. This gives a strict 
quantitative test for deciding whether a particular eigen- 
value reflects a real correlation signal present in the data, 
or is just a spurious noise effect caused by the finite length 
T of the data set. In principle, any eigenvalue signif- 
icantly above the upper cutoff A+ should be treated as 
signal. The variance a of 5xi {t) can be renormalized from 
its starting value cr = 1 by the presence of correlations 
in the data. Indeed, we obtain the best fit of the noise- 
band part of the spectrum for = 0.67, consistent with 
the empirically observed correlations. The three largest 
eigenvalues 2.2,3.5, and 8.7 sufhciently exceed the theo- 
retical upper limit A+ = 1.97 to be attributed to true cor- 
relation patterns. Indeed, in a control test we found that 
the probability that the largest eigenvalue generated by 
an uncorrelated univariate Gaussian signal of the same 
dimensions as our data to exceed 2.2 is around 0.05%. 
This should be contrasted with a typical 5% to 1% con- 
fidence level of correlations between a pair of individual 
indices reported e.g. in In order to check the repro- 
ducibility of largest eigenvalues and their corresponding 
eigenvectors we divided ours set into two consecutive 113- 
point subsets and repeated our analysis. The existence of 
3 outlier eigenvalues did not change, however their val- 
ues have slightly changed. The largest eigenvalue was 
measured to be 9.9 during the first time interval, and 7.8 



during the second. As can be concluded from the inset 
to Fig. 1, and Fig. 2, the corresponding eigenvectors are 
remarkably stable with overlaps between eigenvectors for 
the first and the second subintervals being 0.95, 0.81, and 
0.67 for the largest, the second, and the third eigenvalues 
correspondingly. The largest possible overlap, realized 
when two eigenvectors are identical, is equal to 1. On 
the other hand, overlaps between eigenvectors from the 
noise-band between A_ and A+ seem to be purely random 
(see inset to Fig. 1). Similar effects but with larger num- 
ber of outliers above A+ (up to around 25) were observed 
Q for individual stocks traded on US stock exchanges. 

Components of the three highest ranking eigenvectors, 
measured for our data both in its entirety, and when di- 
vided in two equal subintervals, are given in Table 1 and 
plotted in Fig. 2. The first interesting result is that virtu- 
ally all components of the largest eigenvector are positive, 
which means that there are no indices which are anti- 
correlated with others. Since eigenvectors corresponding 
to different eigenvalues have to be orthogonal to each 
other, other eigenvectors must contain negative compo- 
nents. The first (largest) eigenvector has strong support 
in European and American sectors, while its components 
in the Asian sector are somewhat smaller (yet still pos- 
itive). The second eigenvector, on the other hand, is 
largely dominated by Asian stocks while the third one 
by American stocks. 

Another interesting observation is that all three eigen- 
vector components for some of the Asian emerging mar- 
kets such as China, India, Pakistan, and Taiwan are 
too small to be detected. That means that in the first 
approximation these indices are not influenced by the 
world index dynamics at all. In Europe we saw no such 
correlation-free countries. However, the eigenvector com- 
ponents of Greece, Portugal, Russia, and Turkey were 
somewhat smaller than those of other European stock 
indices. North and South american stock indices have ap- 
proximately equal components with, perhaps, only Peru 
and Venezuela somewhat falling behind. 

It is interesting to compare our findings to those 
which were previously obtained for weeekly returns us- 
ing the Minimal Spanning Tree (MST) technique [||. In 
this work it was also observed that indices are strongly 
grouped by the region with most indices of, say, Asian 
stock indices forming a separate branch of the MST. The 
market indices of Turkey, Greece, India, and Pakistan 
were found to be weakly correlated with the world in- 
dex in both Ref. |^ and our study. We believe that the 
spectral analysis method, used in our work, gives a com- 
plimentary picture of reproducible correlations contained 
in the matrix Cij to that of the MST method. One of the 
strong points about the spectral analysis method is that 
it gives a clear quantitative criterion for separating the 
signal from the noise. Also with just a few large eigen- 
values and their corresponding eigenvectors the output 
of the spectral analysis is easier to interpret and com- 
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pare between, say, different time windows than that of 
the MST technique. 

The leading eigenvector component of a stock market 
index of a given country can serve as a rough measure 
of the level of globalization of the financial sector of this 
country. This point of view is illustrated in Fig. 3 where 
the largest eigenvalue component is plotted as a function 
of the rank of the country in the A.T. Kearney/Foreign 
Policy Magazine Globalization Index"""^ One can 

see a clear correlation between high globalization rank 
(1 being the highest and 50 the smallest) and the leading 
eigenvector component. 

We further decided to explore how the outcome of the 
above eigenvector/eigenvalue analysis dependent on the 
time horizon At over which one computes the returns of 
an index. When instead of daily open-to-close returns we 
repeated our analysis for one day close-to-close returns 
the largest eigenvalues have changed from 2.2, 3.5, and 
8.7 to 2.1, 3.7, and 11.0. A noticeable 25% increase in the 
largest eigenvalue perhaps acn be attributed to markets 
having more time to respond to the news. Also, while 
daily open-to-close returns in Asia, Europe, and Ameri- 
cas have almost no overlap (i.e. time when two markets 
are simultaneously open) , this situation is improved when 
one considers daily close-to-close fluctuations. As shown 
in Fig. 4, the largest eigenvalue continued to grow (al- 
beit slowly) , as the time horizon of close-to-close returns 
was changed from one to ten trading days (weekly re- 
turn usually corresponds to just 5 trading days), reaching 
the value of 16.2 for the longest time horizon. However, 
the largest eigenvectors computed for these very different 
time horizons remained remarkably stable. For exam- 
ple the average overlap between the highest rank eigen- 
vectors computed for these ten different time horizons 
turned out to be 0.99, i.e. these eigenvectors on aver- 
age are only 1% different from each other! Overlaps were 
somewhat smaller for lower ranking eigenvectors with av- 
erage values of 0.77 for the second and 0.61 for the third 
largest eigenvalues. Still, as can be seen in the inset to 
Fig. 4, even in the third eigenvector many of the main 
features are very robust with respect to changes in the 
time horizon. 

In an attempt to establish how relevant are magnitudes 
(as opposed to signs) of price fluctuations to the observed 
correlation patterns we have repeated the above analysis 
using SX,{t) = sign(S'.,(t + At) - Si{t)). The observed 
eigenvectors for different time horizons had 0.99 average 
overlap with those computed using SXi{t) — In Si{t + 
At) — In Si{t). The largest eigenvalue again grew with the 
time horizon from 7.4 for signs of one day close-to-close 
returns to 10.4 for signs of 5-day close-to-close returns. 
Lower rank eigenvectors of the correlation matrix of signs 
also had substantial overlaps of 0.93 and 0.86 with normal 
ones. This allows us to conclude that it is signs and not 
magnitude of returns which are mostly relevant for the 
observed correlation patterns. 



In what was described above we always computed 
(nearly) synchronous correlations of different market re- 
turns on the same day (or the same week for longer time 
horizons). One has to take into account that due to 
the time-zone difference daily open-to-close returns com- 
puted on the same trading day (say, February 13) are 
not actually synchronous, with Asian stock markets in 
the lead, followed by European and later with a small 
overlap by American markets. However, the significance 
of this time zone difference for say weekly returns is much 
less pronounced. To check if the predictability of daily 
returns is restricted to the same trading day or survives 
for several days we have investigated the correlations in 
daily close-to-close returns with all Asian indices shifted 
by S days. The negative values of the shift S corresponds 
to correlations of daily returns of Asian indices \S\ days 
after an observed pattern of returns of the European and 
American stocks, while positive S corresponds to Asian 
indices preceding the rest of the world by S days. The 
simplest way to detect the presence of correlations in this 
case is by calculating the average correlation coefficient 
connecting any of the 14 Asian indices with 15 -I- 8 = 23 
European and American ones. The size of the sample 
over which this average is taken is 14 x 23 = 322. To 
check for reproducibility of the observed patterns we di- 
vided our data set consisting of 226 trading days into 
three equal length segments and calculated this average 
for each segment independently. The results are shown 
in the main panel of Fig. 5. Reproducible correlations 
seem to survive for up to 3 days on the negative part of 
the axis, corresponding to the reaction of Asian stocks to 
changes in European and American indices. At least the 
sign of the average correlation coefficient was found to be 
consistently positive in all three of our segments. 5 = — 1 
has the largest magnitude of correlations. These (posi- 
tive) correlations represent the next trading day reaction 
of Asian stocks to changes in European and American 
ones. Due to time-zone differences a somewhat smaller 
(yet still positive) correlation coefficients observed for 
S — correspond to the opposite effect, i.e the response 
of European and American stocks to the events in Asia 
on the day before. On the positive S side there seems to 
be a reproducible negative correlation at S = 1 followed 
by a noisy signal for larger values of S. Our results for 
S — and S — 1 are in agreement with positive corre- 
lations between open-to-close returns at Tokyo and New 
York stock exchanges that were previously reported in 
the economics literature Q as well as in recent econo- 
physics papers The inset to Fig. 5 show the results 
of the same analysis for the whole 226-day interval when 
either Asian, or European and Asian stocks are shifted 
by S days. It is interesting that in this graph which has 3 
times better statistics than the main panel of the figure, 
there seem to be oscillations of correlation coefficients 
with period of two trading days for negative values of S 
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FIG. 1. The histogram of eigenvalues of the matrix dj. The straight line is a fit with the Eq. (1) using T — 226, — 37, 
and various values of a, corresponding to the exclusion of fluctuations contained in the largest (cr^ = 0.76), or two largest 
{a^ = 0.67) eigenvalues. The inset shows the overlap of eigenvectors computed for two consecutive 113-day intervals as a 
function of the rank of an eigenvalue. The three (perhaps even four) leading eigenvectors clearly have a higher than random 
overlap. 
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FIG. 3. The component of the highest ranking eigenvector as a function of the globaUzation rank of the country from Ref.[l]. 
The straight hne is a hnear fit to the data. 
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FIG. 4. Three largest eigenvalues of the correlation matrix of close-to-close returns as a function of the time horizon (number 
on days used to calculate the returns). The inset shows the components of the third eigenvector for all ten time horizons. 
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FIG. 5. The mean value of the correlation coefficient connecting Asian indices to the rest of the world as a function of the 
shift S. Negative values of S correspond to Asian indices taken \S\ days later than the rest of the world. Three data sets were 
taking in three equal length subintervals of our data set. The inset shows the same analysis repeated including all data points 
with Asian indices (circles) and Asian and European indices (squares) shifted. Note oscillations for negative S. 
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Country 


Index name, Symbol 








2nd, e.v. 


2nd, e.v. 




3rd, e.v. 


3rd, e.v. 


3rd, e.v. 


Rank 








all days 


1 St half 


2nd half 


all days 


1st half 


2nd half 


all days 


1 St half 


2nd half 


in [1] 


1 


Austria 


ATX, ■^ATX 


0.20 


0.22 


0.12 


0.02 


0.03 


0.03 


-0.19 


-0.19 


-0.30 


7 


2 


BslQium 


BEL-20 , '*BFX 


0.16 


0.19 


0.1 1 


-0.17 


-0.21 


-0.1 1 


-0.12 


-0.19 


-0.04 


N/A 


3 


DGnmark 


KFX , '*KFX 


0.21 


0.21 


0.21 


0.07 


0.06 


0.1 1 


-0.17 


-0.1 1 


-0.23 


1 1 


4 


Finland 


HELSINKI GENERAL , "HEX 


0.22 


0.25 


0.25 


0.02 


-0.05 


0.07 


-0.14 


-0.08 


-0.13 


5 


5 


France 


CAC40 , ''FCHI 


0.21 


0.21 


0.22 


-0.28 


-0.28 


-0.26 


-0.15 


-0.06 


-0.19 


16 


6 


GGrmany 


DAX , "GDAXI 


0.27 


0.25 


0.28 


-0.14 


-0.15 


-0.1 1 


-0.14 


-0.09 


-0.17 


14 


7 


Grssce 


GENERAL SHARE , MTG 


0.10 


0.15 


-0.02 


0.07 


0.07 


0.04 


0.10 


0.00 


0.16 


24 


8 


Nstherlands 


AEX GENERAL , '^AEX 


0.23 


0.22 


0.22 


-0.25 


-0.24 


-0.25 


-0.17 


-0.14 


-0.19 


2 


9 


Norway 


TOTAL SHARE , '^NTOT 


0.22 


0.20 


0.25 


0.07 


0.03 


0.14 


-0.12 


-0.1 1 


-0.15 


9 


10 


Portugal 


BVL30 , '^BVLSO 


0.14 


0.12 


0.17 


-0.26 


-0.29 


-0.21 


-0.07 


-0.07 


-0.01 


15 


11 


Russia 


MOSCOW TIMES , '^MTMS 


0.12 


0.10 


0.15 


0.05 


0.1 1 


-0.06 


-0.07 


-0.14 


0.01 


45 


12 


Sweden 


STOCKHOLM GENERAL , ''SFOG 


0.26 


0.25 


0.27 


0.01 


0.02 


0.01 


-0.17 


-0.13 


-0.20 


3 


13 


Switzerland 


SWISS MARKET , "SSMt 


0.21 


0.24 


0.12 


-0.23 


-0.23 


-0.19 


-0.03 


-0.01 


-0.16 


4 


14 


Turkey 


ISE NATIONAL-100 , '^XUIOO 


0.1 1 


0.1 2 


0.08 


0.17 


0.16 


0.17 


-0.08 


-0.17 


0.12 


37 


15 


UK 


FTSE100 . 'FTSE 


0.26 


0.25 


0.29 


0.01 


0.05 


0.06 


0. 17 


0. 15 


0.1 1 


8 


16 


Australia 


ALL ORDINARIES . "'AORD 


0. 1 1 


0.1 1 


0. 1 1 


0.35 


0.28 


0.39 


-0.02 


-0.06 


0.07 


23 


"17 




SHANGHAI COMPOSITE , "SSEC 


0.00 


0.01 


-0.03 


0.07 


0.06 


0.08 


-0.04 


-0.01 


-0.12 


48 


18 


Honk Kong 


HANG SENG , '^HSE 


0.13 


0.1 1 


0.16 


0.25 


0.23 


0.23 


-0.17 


-0.19 


-0.03 


N/A 


19 




BSE30 , ''BSESN 


0.02 


0.02 


0.03 


0.07 


0.21 


-0.07 


0.03 


-0.07 


0.12 


49 


20 


Indonesia 


JAKARTA COMPOSITE , '^JKSE 


0.06 


0.07 


0.03 


0.24 


0.28 


0.18 


0.05 


0.02 


0.04 


38 


21 


Japan 


NIKKEI 225 , '^N225 


0.08 


0.07 


0.10 


0.30 


0.22 


0.36 


-0.07 


-0.07 


-0.01 


29 


22 


Malaysia 


KLSE COMPOSITE , '^KLSE 


0.07 


0.06 


0.09 


0.15 


0.19 


0.01 


-0.17 


-0.28 


0.08 


20 


23 


New Zealand 


NZSE40 , '^NZ40 


0.07 


0.09 


0.04 


0.29 


0.25 


0.32 


0.00 


0.01 


-0.01 


21 


24 


Pakistan 


KARACHI100,''KSE 


0.00 


0.00 


-0.01 


0.00 


0.03 


-0.02 


0.08 


0.08 


-0.08 


N/A 


25 


Philippines 


PSE COMPOSITE , -^PSI 


0.07 


0.06 


0.09 


0.17 


0.19 


0.13 


0.18 


0.27 


-0.01 


34 


26 


Singapore 


STRAITS TIMES , "STI 


0.13 


0.10 


0.19 


0.24 


0.25 


0.19 


0.03 


0.03 


0.11 


1 


27 


South Korea 


SEOUL COMPOSITE , ^KS1 1 


0.09 


0.07 


0.10 


0.19 


0.15 


0.23 


-0.03 


0.10 


-0.11 


31 


28 


Taiwan 


TAIWAN WEIGHTED , ^TWII 


0.00 


0.01 


0.00 


-0.03 


-0.01 


-0.08 


-0.09 


-0.07 


-0.08 


N/A 


29 


Thailand 


SET , ^SETI 


0.12 


0.11 


0.12 


0.08 


0.14 


-0.03 


-0.07 


-0.06 


-0.02 


30 


30 


Argentina 


MERVAL.^MERV 


0.21 


0.21 


0.19 


0.01 


0.08 


-0.08 


0.31 


0.27 


0.31 


39 


31 


Brazil 


BOVESPA , ^BVSP 


0.19 


0.17 


0.20 


-0.04 


-0.01 


-0.08 


0.42 


0.44 


0.32 


44 


32 


Canada 


TSE 300 COMPOSITE , ^TSE 


0.23 


0.25 


0.20 


-0.06 


0.01 


-0.14 


0.17 


0.14 


0.16 


10 


33 


Chile 


IPSA , /^IPSA 


0.20 


0.19 


0.19 


-0.01 


-0.02 


0.02 


0.30 


0.29 


0.30 


26 


34 


Mexico 


IPC , '^MXX 


0.20 


0.18 


0.22 


-0.06 


-0.08 


-0.03 


0.29 


0.26 


0.30 


41 


35 


Peru 


LIMA GENERAL , '^IGRA 


0.17 


0.19 


0.11 


0.14 


0.17 


0.11 


0.17 


0.10 


0.20 


43 


36 


USA 


SP500 , '^SPC 


0.21 


0.20 


0.21 


-0.13 


-0.07 


-0.21 


0.31 


0.32 


0.26 


12 


37 


Venezuela 


IBC , '^IBC 


0.13 


0.12 


0.12 


0.10 


0.13 


0.07 


0.09 


0.04 


0.04 


35 



TABLE I. Components of three leading eigenvectors computed for the whole 226-day time interval and its first and second 
halves. The last column is the rank of the globalization index of the country as defined in Ref. [1] 
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